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An  ab  Initio  Study  of  the  Kinetics  of  the  Reactions  of  Halomethanes  with 
the  Hydroxyl  Radical.  Part  1 :  CH2Br2 

Florent  LOUIS*,  Robert  E.  HUIE  and  Michael  J.  KURYLO 

Physical  and  Chemical  Properties  Division,  National  Institute  of  Standards  and  Technology, 
Gaithersburg,  Maryland  20899,  USA 

Abstract 

Ab  initio  calculations  using  the  Moller-Plesset  perturbation  theory  were  carried  out  on 
the  H-atom  abstraction  reaction  from  dibromomethane  by  hydroxyl  radical  attack.  Geometry 
optimization  and  vibrational  frequency  calculations  at  the  MP2  level  were  performed  on  all 
reactants,  products  and  the  transition  state  with  both  the  6-3 1  lG(d,p)  and  6-3 1  lG(2d,2p)  basis 
sets.  Utilizing  the  geometry  parameters  optimized  at  the  MP2/6-31  lG(2d,2p)  level  of  theory, 
single-point  energy  calculations  were  carried  out  with  increasing  basis  set  sizes,  from  6- 
31  lG(2d,2p)  to  6-31 1++G(3df,3pd)  and  with  both  the  MP2  and  MP4SDTQ  methods. 
Canonical  transition  state  theory  was  used  to  predict  the  rate  constants  as  function  of  the 
temperature  (240-380  K).  For  the  kinetic  parameters  of  this  reaction,  reasonable  agreement 
with  the  experimental  values  was  obtained  at  the  PMP2/6-3 1 1  G(3df,2p)//MP2/6-3 1 1  G(2d,2p) 
level  of  theory.  With  the  largest  basis  set,  6-3 1 1++G(3df,3pd),  the  calculated  rate  constants 
are  in  very  good  agreement  with  their  experimental  counterparts. 

*  Corresponding  author,  Fax  :  (301)  975-3672;  E-mail:  flouis@nist.gov 
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Introduction 

In  this  paper,  we  initiate  a  study  of  the  application  of  ab  initio  transition  state  theory  to 
the  kinetics  of  hydrogen-atom  abstraction  by  the  hydroxyl  radical.  The  objective  of  this  study 
is  to  calculate  the  enthalpy  of  the  reaction  for  molecules  containing  various  functional  groups, 
most  importantly  those  containing  halogen  atoms,  and  to  calculate  the  rate  constants  as  a 
function  of  temperature.  From  the  latter,  we  calculate  Arrhenius  A-factors  and  activation 
energies.  In  the  present  paper,  we  report  the  results  of  our  investigation  of  the  reaction: 

OH  +  CH2Br2  -»  H20  +  CHBr2  (1) 

The  initial  objective  is  to  establish  that  level  of  theory  (both  method  and  basis  set)  which  is 
sufficient  to  describe  the  kinetics  for  this  reaction.  The  computational  difficulties  in  treating 
the  large  electronic  system  associated  with  bromine  make  this  a  serious  test  of  this  use  of 
theory  and  the  computational  level  optimum  for  this  reaction  should  be  adequate  for  other 
molecules  containing  bromine  or  lighter  halogens.  Further,  bromine  is  an  active  agent  in  fire 
suppression  and  inhibition.1  In  the  screening  for  suitability  as  a  fire  suppressant,  the 
reactivity  of  the  agent  is  an  important  consideration,  both  with  regards  to  its  utility  and  its 
environmental  acceptability. 

There  have  been  some  early  applications  of  transition-state  theory  to  the  reactions  of 
OH  with  haloalkanes,  making  use  of  thermochemical  kinetic  tools,  which  provided  estimates 
of  the  Arrhenius  pre-exponential  factors.2*3  A  few  ab  initio  studies  also  have  been  carried  out 
on  the  reactions  of  hydroxyl  radicals  with  halogen-substituted  methanes4'6  or  ethanes.2'1 1 
None  of  these,  however,  have  been  for  bromine-substituted  methanes.  In  addition,  there  have 
been  several  ab  initio  studies  on  the  reactions  of  fluorine  and  chlorine  atoms  with  halogenated 
methanes12'14  and  ethanes13*16  plus  one  on  methyl  radical  reactions  with  halomethanes17 
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and  a  theoretical  study  of  the  enthalpies  of  formation  of  fluoromethanes.18  Recently,  two  ab 
initio  studies  dealing  with  the  structures,  vibrational  frequencies,  thermodynamic  properties, 
and  bond  dissociation  energies  of  bromomethanes  and  bromomethyl  radicals  were 

published.  19-21 

Table  I  summarizes  the  available  results  of  the  experimental  studies  of  the  reaction  of 
OH  with  CH2Br2.  The  rate  constant  was  measured  using  absolute  techniques  by  Mellouki  et 
al.22  and  Zhang  et  al.22  as  well  as  by  relative  techniques  by  DeMore24  and  Orlando  et  al.25  it 
was  also  measured  recently  in  our  laboratory  at  298  K.2^  with  the  exception  of  the  relative 
results  of  DeMore  24,  which  have  been  recalculated  with  the  recommended  value  for  the 
reference  reaction  of  OH  with  CH2C12,27  the  measurements  all  agree  with  the  value  k(298)  = 
1.2  x  10-13  cm3  molecule-1  s-1. 

The  kinetic  parameters  for  the  reaction  of  OH  with  CH2Br2  are  consistent  with  a 
reaction  that  takes  place  by  hydrogen  abstraction.  Although  recent  work  has  demonstrated  that 
the  reaction  of  the  hydroxyl  radical  with  CF3I  proceeds  rapidly  by  I-atom  abstraction.28  The 
lack  of  reactivity  of  CF3Br  towards  OH27  argues  against  any  significant  contribution  of  a 
bromine  abstraction  reaction  pathway  for  OH  +  CH2Br2. 

Computational  Methods29 

All  calculations  described  below  were  carried  out  with  the  program  package  Gaussian 
94^0  on  Cray  C90/6256  and  SGI  Origin  2000  computers.  Fully  optimized  geometry 
parameters  and  corresponding  vibrational  frequencies  for  the  reactants,  products  and  the 
transition  state  were  computed  with  the  second  order  Moller-Plesset  perturbation  theory 
(UMP2)  with  the  6-3 1 1  G(d,p)  and  6-3 1 1  G(2d,2p)  basis  sets.  Spin  contamination  in  the  open- 
shell  systems  was  taken  into  account  by  a  projection  method  included  in  the  Gaussian  94 
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program  and  is  noted  in  the  results  as  PMPn  (n  =  2  or  4). 3 1.32  Vibrational  frequencies  have 
been  used  for  the  characterization  of  stationary  points,  zero-point  energy  (ZPE)  correction  and 
transition  state  theory  (TST)  computations  of  the  reaction  rate  constant.  The  calculated 
transition  state  was  confirmed  to  have  one  imaginary  frequency.  All  the  energies  quoted  and 
discussed  in  the  present  paper  include  the  zero-point  energy  correction  with  unsealed 
vibrational  frequencies. 

The  optimized  geometry  parameters  at  the  MP2/6-31  lG(2d,2p)  level  of  theory  were 
utilized  in  the  calculations  of  single-point  energies.  These  calculations  were  carried  out  at 
both  the  MP2  and  MP4SDTQ  levels  of  theory  and  with  basis  set  sizes  increasing  from  6- 
31  lG(2d,2p)  to  6-31 1++G(3df,3pd).  Canonical  transition  state  theory  and  tunneling 
corrections  were  used  to  predict  the  rate  constant  over  the  same  range  of  temperatures  as  the 
experimental  measurements  (240-380  K).33 

Results  and  Discussion 

1.  Geometry  Parameters  and  Vibrational  Frequencies 

1.1.  Geometry  parameters 

Reactants  and  products'.  Table  II  shows  optimized  geometry  parameters  at  different 
levels  of  theory  and  also  the  available  experimental  data  characterizing  the  structures  of  these 
species.  All  the  tested  theoretical  models  produce  acceptable  parameters  for  these  species. 

Transition  state :  The  optimized  parameters  for  the  transition  state  at  the  two  levels  of 
•  theory  are  shown  in  Table  III  and  Figure  1.  A  comparison  of  results  obtained  with  the  MP2 
method  and  the  6-31  lG(d,p)  and  6-31  lG(2d,2p)  basis  sets  shows  that  the  parameters 
calculated  for  the  geometry  of  the  transition  state  structure  are  very  similar.  These  results  are 
reflected  in  the  transition  state  parameter  L  which  is  defined  as  the  ratio  between  the  increase 
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in  length  of  the  C-H  bond  being  broken  and  the  elongation  of  the  O-H  bond  being  formed, 
each  with  respect  to  its  equilibrium  value  in  the  reactant  (CH2Br2)  and  the  product  (H20)J2 
This  parameter  characterizes  the  most  important  aspect  of  the  geometric  structure  of  the 
transition  state.  When  the  basis  set  included  two  d-orbital  polarization  functions  on  the  heavy 
atoms  (C  and  Br)  and  two  p-orbital  polarization  functions  on  the  H  atoms,  the  parameter  L 
changes  only  slightly.  This  suggests  that  the  MP2/6-31  lG(2d,2p)  level  of  theory  gives 
adequate  results  for  the  geometry  of  the  transition  state.  Consequently,  these  MP2/6- 
31  lG(2d,2p)  optimized  geometry  parameters  will  be  used  for  the  single-point  energy 
calculations. 

1.2.  Vibrational  frequencies 

The  vibrational  frequencies,  unsealed  zero-point  energies  and  internal  energy 
corrections  at  298  K  for  reactants,  transition  state  and  products  are  given  in  Tables  IV  and  V. 
For  reactants  and  products,  calculated  vibrational  frequencies  values  at  the  MP2/6-3 1  lG(d,p) 
and  MP2/6-31  lG(2d,2p)  are  about  7%  greater  than  the  experimental  values.  In  the  transition 
state,  the  eigenvector  corresponding  to  the  imaginary  frequency  is  primarily  a  motion  of  the 
reactive  hydrogen  atom  being  transferred  between  the  C  and  O  centers.  The  calculated  value 
of  the  imaginary  frequency  is  about  6%  larger  at  the  MP2/6-31  lG(2d,2p)  level  of  theory  than 
at  the  MP2/6-3 1  lG(d,p)  level.  This  fact  was  previously  observed  for  the  series  of  H-atom 
abstraction  reactions  F  +  CHC13.XFX  (x  =  0, 1, 2  or  3).^ 

2.  Reaction  enthalpies  and  activation  energies 

The  effect  of  the  basis  set  and  the  influence  of  higher-order  electron  correlation  and 
spin  projections  on  the  calculated  total  energies  of  the  transition  state,  reactants  and  products 
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are  shown  in  Table  VI.  All  the  single-point  energy  calculations  were  done  with  the  MP2/6- 
31  lG(2d,2p)  optimized  geometry  parameters  except  for  the  first,  which  used  the  MP2/6- 
31  lG(d,p)  optimized  geometry  parameters. 

In  the  first  step,  we  calculated  single-point  energies  for  reactants,  products  and 
transition  state  at  different  levels  of  theory.  Corresponding  reaction  enthalpies  (AH)  and 
activation  energies  (Ea),  taking  into  account  the  zero-point  energy  difference  and  thermal 
energy  corrections  at  298  K  (calculated  using  the  vibrational  frequencies  given  in  Table  IV), 
are  displayed  in  Table  VII  (columns  2  and  5). 

The  two  geometry  optimization  levels  of  theory  give  calculated  reaction  enthalpies  that 
were  very  close  to  the  literature  values  based  on  A/H°  at  298  K  for  OH,^4  CH2Br2,35  H20,^4 

and  CHBr2,36  particularly  when  the  stated  uncertainties  for  CH2Br2  (±9  kJ  mol'1  )  and  CHBr2 
(+5  kJ  mol"1)  are  taken  into  consideration.  Adding  diffuse  functions  to  the  basis  set  with  MP2 
method  led  to  calculated  reaction  enthalpies  AH  outside  the  experimental  uncertainties.  With 
all  the  other  tested  theoretical  models,  a  satisfying  agreement  was  obtained  between 
experimental  reaction  enthalpies  and  their  calculated  counterparts. 

In  order  to  get  a  much  closer  agreement  between  experimental  and  ab  initio  values, 
especially  for  the  basis  sets  including  diffuse  functions,  reaction  enthalpies  for  reaction  (1) 
were  also  corrected  for  the  different  levels  of  theory  by  utilizing  the  isodesmic  reaction 

CH3  +  CH2Br2  -»  CH4  +  CHBr2  (2) 

to  calculate  D298  (H-CHBr2).  The  direct  calculation  of  this  bond  dissociation  energy  from  the 
reaction  CH2Br2  — >  H  +  CHBr2  is  expected  to  be  systematically  in  error  due  to  insufficient 
treatment  of  electron  correlation  and  incompleteness  of  the  basis  sets.  The  use  of  an 
isodesmic  reaction,  such  as  reaction  (2),  provides  an  indirect  method  which  may  lead  to  a 
more  reliable  value  of  the  bond  dissociation  energy.  The  specific  property  of  an  isodesmic 
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reaction  is  that  there  are  the  same  number  and  type  of  bonds  on  each  side  of  the  equation,  so 
that  the  errors  mentioned  above  are  largely  canceled  in  the  calculated  value  of  D298  (H- 
CHBr2).  The  reaction  enthalpy  for  reaction  (2)  was  calculated  at  the  same  levels  of  theory  as 
that  for  reaction  (1).  Taking  the  experimental  value  of  D298(H-CH3)  =  440  kJ  mol'1, 3 7  we 
obtain  calculated  D298  (H-CHBr2)  values  at  the  various  levels  of  theory  (Table  VII,  column  3). 
Then  the  reaction  enthalpy  for  reaction  (1)  are  determined  as  A,!!  (ISO)  =  D298(H-CHBr2)  - 
D298  (H-OH).  Table  VII  (column  4)  contains  these  new  values  which  can  be  compared  to  the 
direct  calculations  (column  2)  and  to  their  experimental  counterparts.  The  improvement 
observed  is  quite  significant,  especially  for  the  basis  sets  with  diffuse  functions. 

The  activation  energy  calculated  at  the  PMP2/6-31  lG(d,p)//MP2/6-311G(d,p)  level  of 
theory  is  about  13-14  kJ  mol'1  higher  than  the  measured  value.  Improvement  of  the  basis  set  to 
6-31  lG(2d,2p)  yields  a  decrease  in  the  value  activation  energy  by  2.9  kJ  mol'1.  Increasing  the 
basis  set  size  by  adding  more  polarization  functions  and  diffuse  functions  led  to  decrease  the 
values  of  activation  energy  so  that  the  differences  between  calculated  and  experimental 
activation  energies  for  the  OH  +  CH2Br2  reaction  become  smaller,  i.e.  5.0  to  3.5  kJ  mol"1 
respectively  for  the  PMP2/6-31  lG(3df,2p)  and  the  largest  basis  set  with  both  MP2  and 
MP4STDQ  methods. 

3.  Kinetic  parameters  and  tunneling  effect 

The  knowledge  of  the  geometry  and  vibrational  frequencies  characterizing  the 
transition  state  (Tables  III  and  IV)  enables  us  to  calculate  the  corresponding  activation  entropy 
(A*S)  and  the  Arrhenius  A- factor  using  conventional  transition  state  theory.  Table  VIII 
collects  the  values  of  molar  entropies  calculated  at  298  K  for  reactants  and  the  transition  state, 
the  ones  of  resulting  activation  entropies  and  the  Arrhenius  A-factors  calculated  at  the  two 
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geometry  optimization  levels.  At  both  the  MP2/6-31  lG(d,p)  and  MP2/6-31  lG(2d,2p)  levels 
of  theory,  good  agreement  with  the  calculated  and  experimental  values  of  the  molar  entropies 
were  obtained.  The  electronic  contribution  of  the  molar  entropy  of  the  hydroxyl  radical  took 
into  account  the  multiplicity  of  the  states  and  the  energy  gap  of  139.7  cm*1  between  the  2Ylm 
and  2Tll/2  levels. 34  The  resulting  A-factors  are  in  good  agreement  with  experiment. 

For  the  reactions  involving  a  hydrogen-atom  transfer,  the  tunneling  effect  must  be 
estimated.  The  great  importance  of  this  correction  suggests  that  a  preferable  approach  to 
obtaining  the  Arrhenius  parameters  is  through  the  calculation  of  the  rate  constants  for  the 
reaction  over  a  range  of  temperatures,  with  a  subsequent  fit  of  the  rate  constants  to  the 
Arrhenius  equation.  The  rate  constants  were  calculated  over  the  same  temperature  range  as  the 
measured  values  (i.e.  240  to  380  K)  using  the  canonical  transition  state  theory  formalism^ 
including  the  symmetrical  Eckart  tunneling  correction,  r(T)38 


kbT 

k(T)=T(T)  x  —7—  x  — 

n  Q 


Q*(T) 

OH  CO  QcH2Br2 


where  Q*(T),  QOH(T)  and  are  the  total  partition  functions  for  the  transition  state, 

hydroxyl  radical,  and  dibromomethane  at  temperature  T;  AE  is  the  activation  energy  including 
thermal  corrections  to  the  internal  energy  and  zero-point  energy;  kb  is  the  Boltzman  constant; 
and  h  the  Planck  constant.  Calculated  tunneling  parameters  using  the  symmetrical  Eckart’s 
correction  are  included  in  the  range  4.0  to  10.0  as  function  of  the  temperature.  At  room 
temperature,  their  values  are  about  5.8  to  7.2  depending  on  the  level  of  theory  used. 

Symmetrical  Eckart  tunneling  is  more  appropriate  to  calculate  the  tunneling  factor  than 
the  Wigner’s  correction. 38  For  even  better  agreement,  more  sophisticated  methods  such  as 
the  small  curvature  approximation  are  necessary. 3 9  in  this  case,  detailed  knowledge  of  the 


9 


reaction  path  is  needed  and  this  method  is  very  time  consuming  with  a  large  electronic  system 
such  as  that  for  the  bromine  atom. 

The  results  of  the  calculations  of  the  rate  constants  at  different  levels  of  theory  are 
shown  in  Figure  2  and  a  summary  of  the  calculated  Arrhenius  parameters  including  the 
tunneling  correction  is  given  in  Table  IX.  Regardless  of  the  level  of  theory,  the  calculated 
Arrhenius  A-factors  are  similar  and  in  reasonable  agreement  with  the  experimental  value.  The 
calculated  values  of  the  activation  energy,  however,  depend  strongly  on  the  level  of  theory. 
With  the  largest  basis  set,  6-31 1++G(3df,3pd),  and  with  either  the  MP2  or  MP4SDTQ 
methods,  the  predicted  rate  constants  are  in  very  good  agreement  with  their  experimental 
counterparts  over  the  whole  temperature  range.  At  the  PMP2/6-3 1  lG(3df,2p)//MP2/6- 
31  lG(2d,2p)  level  of  theory,  which  is  significantly  less  computationally  demanding,  the 
calculated  rate  constants  are  47  %  lower  than  the  measured  values  over  the  temperature  range 
(240-380)  K. 


Conclusion 

Ab  initio  calculations  were  performed  at  different  levels  of  theory  for  the  H-atom 
abstraction  reaction  between  dibromomethane  and  hydroxyl  radicals.  The  geometry 
parameters  for  reactants,  products  and  transition  state  were  first  fully  optimized  at  the  MP2/6- 
3 1  lG(d,p)  level  of  theory  and  then  reoptimized  with  the  6-3 1  lG(2d,2p)  basis  set.  The 
parameter  L  changed  only  slightly,  suggesting  that  the  MP2/6-3 1  lG(2d,2p)  level  of  theory 
gives  adequate  results  for  the  geometry  of  the  transition  state.  Hence  it  appears  that  no  higher 
basis  set  is  needed  for  the  optimization  of  the  geometry  parameters  in  subsequent 
investigations  of  the  reactions  of  OH  with  other  halomethanes.  The  calculation  of  the 
energetics  of  the  reaction  and,  subsequently,  the  kinetic  parameters,  proved  more  dependent 
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on  the  level  of  theory,  particularly  on  the  nature  and  extent  of  the  basis  set.  The  best  fit  to  the 
experimental  data  was  with  the  6-31 1++G(3df,3pd)  basis  sets,  employing  either  the  MP2  or 
MP4SDTQ  methods.  A  quite  reasonable  agreement,  within  a  factor  of  two,  was  obtained  with 
the  far  less  computationally  intensive  MP2  method  with  a  6-31  lG(3df,2p)  basis  set.  These 
levels  of  theory  are  applied  in  our  subsequent  work  to  a  broader  range  of  halocarbon  reactions 
with  the  hydroxyl  radical.  The  results  for  reactions  involving  halomethanes  of  type  CHXYZ 
(X,  Y,  Z  =  H,  F,  Cl  or  Br)  indicate  that  the  energetics  and  kinetics  are  well  described  with 
these  levels  of  theory  and  that  the  calculated  rate  constants  are  in  very  good  agreement  with 
their  experimental  counterparts. 
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Table  I :  Experimental  rate  parameters  for  the  reaction  OH  +  CH2Br2  -»  H20  +  CHBr2 


A 

(cm3  .molecule'1  .s'1) 

E, 

(kJ  mol'1) 

T  (K) 

k  (298  K) 

(cm3.molecule'1.s1) 

Method 

Ref. 

1.9  x  10'12 

7.0 

244-370 

l.l4x  10'13 

FP/LIFa 

22 

3.2  x  10'12 

7.5 

293-375 

1.5  x  10'13 

RRb 

24 

298 

1.2  x  10'13 

RRb 

25 

1.5  x  10'12 

6.0 

288-368 

1.3  x  10'13 

DF/RFC 

23 

a  tt'tv  rr  rr  .  t-*i _ 1_  Til _ 

_ i  —  rr _ - 

298 

-  T _ 1 _ 1  T71 

1.2  x  lO-13 

i _  b  nn  .  r> , 

FP/RFd 

.1  Ti  ~  a  ~ 

26 

C  T\T?/T»T?  . 

a  FP/LIF  :  Flash  Photolysis/Laser  Induced  Fluorescence,  b  RR  :  Relative  Rate,  c  DF/RF  : 
Discharge  Flow/Resonance  Fluorescence, 4  FP/RF  :  Flash  Photolysis/Resonance  Fluorescence 
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Table  II  :  Optimized  geometry  parameters  for  reactants  and  products  involved  in  H- 


abstraction  reaction  of  CH2Br2  with  hydroxyl  radicals 


Parameter* 

MP2 

6-311G(d,p) 

MP2 

6-311G(2d,2p) 

Expt 

OH 

r  (OH) 

0.967 

0.965 

0.97  lb 

CH2Br2 

r  (CH) 

1.086 

1.078 

1.097c 

r  (CBr) 

1.929 

1.933 

1.925 

0  (HCH) 

111.6 

111.9 

110.9 

0  (BrCH) 

107.8 

107.9 

<j>  (BrCHH) 

118.2 

118.5 

h2o 

r  (OH) 

0.958 

0.957 

0. 958d 

0  (HOH) 

102.4 

103.3 

104.5 

CHBr2 

r  (CH) 

1.082 

1.074 

r  (CBr) 

1.860 

1.861 

0  (BrCH) 

116.0 

116.1 

<j>  (BrCHBr) 

148.4 

149.2 

*  Bond  lengths  are  in  Angstroms,  bond  angles  0  and  dihedral  angles  §  are  in  degrees. 


Experimental  values  taken  from b  ref.  ^4, c  ref.  40 , d  ref.  4! 
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Table  III  :  Optimized  geometry  parameters  for  transition  state  involved  in  H-atom 


abstraction  reaction  of  CH2Br2  with  hydroxyl  radicals 


Parameter* 

MP2 

6-311G(d,p) 

MP2 

6-311G(2d,2p) 

r  (C-Hr) 

1.192 

1.188 

r  (CH) 

1.087 

1.078 

r(CBr(a)) 

1.910 

1.911 

r(CBr(b)) 

1.917 

1.919 

rCOHJ 

1.296 

1.296 

r(HO) 

0.968 

0.967 

0  (HCHr) 

108.2 

108.2 

0  (Br^CH*) 

105.8 

106.0 

0  (Br^CH,) 

108.0 

107.8 

0  (OHrC) 

166.2 

165.3 

0  (HOHJ 

97.1 

97.2 

♦  (Br(a)CHRH) 

117.8 

118.1 

<1)  (B^CHrH) 

-118.3 

-118.4 

4>  (ohrch) 

48.5 

51.9 

<j>  (HOHrC) 

50.7 

51.5 

Lb 

0.314 

0.324 

a  Bond  lengths  are  in  Angstroms,  bond  angles  0  and  dihedral  angles  <|)  are  in  degrees  ;  the 


hydrogen  atom  involved  in  H-atom  abstraction  is  noted  HR,  (a)  refers  to  the  in-plane  bromine 
and  (b)  to  the  out-of-plane  bromine 

b  the  parameter  L  is  the  ratio  between  the  elongation  value  of  the  C-H  bond  and  the  elongation 
of  the  O-H  bond  :  L  =  5r(CH)/8r(OH) 
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Table  IV  :  Calculated0'2*  vibrational  frequencies  (cm'1)  for  reactants,  transition  state  and 
products  involved  in  H-atom  abstraction  reaction  of  CHzBrz  with  hydroxyl  radicals 

Species  Vibrational  frequencies  (cm1) 

__  0X3353 

(2)  3833 
3735a 

CHzBrz  (1)  178,  604,  692,  ,835,  1150, 1265, 1459,  3169,  3254 
(2)  175,  590,  672,  833,  1141,  1239,  1458,  3183,  3270 
169,  588,  653,  812,  1095,  1195,  1382,  3009,  3073b 
TS  (1)  2075i,  70, 102, 155,  183,473,651,713,805,  959,  1210, 1319, 1405,  3211,3840 

(2)  2200i,  76, 109, 156,  182, 463,  640,  698,  803,  957, 1 191, 1307, 1395,  3225,  381 1 

H^O  (1)  1667, 3905, 4013 

(2)  1685,3875,  3989 
1595,  3657,  3756b 

CHBrz  (1)  193,  492,  653,  805,  1241,  3267 
(2)  189,  486,  644,  790,  1214,  3279 
633,  778,  1164° 

(,)  MP2/6-31  lG(d,p)  and (2)  MP2/6-31  lG(2d,2p) 

The  experimental  values  of  the  vibrational  frequencies  are  in  italics.  These  values  are  taken 
from  a  ref  $4  b  ref.  42  c  ref  43 
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Table  V  :  Calculated  zero-point  Energy  and  thermal  energy  corrections  at  298  K  for 
reactants,  transition  state  and  products  involved  in  H-atom  abstraction  reaction  of 
CH2Br2  with  hydroxyl  radicals 


Zero-point  Energy* 

Thermal  Energy* 

OH 

23.1b/22.9c 

6.2 

CH2Br2 

75.4/75.1 

10.1 

H20 

57.3/57.1 

7.5 

CHBr2 

39.8/39.5 

10.2 

TS 

90.3/89.8 

16.4 

a  Units  are  kJ  mol'1;  geometry  optimization  calculations  at  bMP2/6-311G(d,p)  and  cMP2/6- 


31  lG(2d,2p)  levels  of  theory 


Table  VI :  Total  energies  (Hartree)  for  species  involved  in  the  OH  +  CH2Br2  reaction 
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All  the  geometry  optimization  calculations  were  done  with  the  MP2/6-3 1  lG(2d,2p)  except  the  first  with  MP2/6-3 1  lG(d,p) 


Table  VII  :  Reaction  enthalpies  A.H,  bond  strengths  D298(H-CHBr2),  reactions  enthalpies  A,!!  (ISO)  calculated  using  an  isodesmic 


of  thermal  energies  (AZPE  +  thermal  energy  corrections), c  including  the  sum  of  thermal  energies  (AZPE  +  thermal  energy  corrections  +  RT), d  in 
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Table  VIII :  Entropies  of  reactants  (S),  transition  state  (S')  and  activation  entropy  (A'S) 


and  preexponential  factor  calculated  at  298  K  at  different  levels  of  theory 


S 

J  mol 1  K 1 

S' 

J  mol1  K 1 

A'S 

J  mol 1  K 1 

A  x  1012 

cm3  molecule'1  s'1 

MP2/6-311G(d,p) 

MP2/6-31  lG(2d,2p) 

OH 

185.8 

185.8 

(183.7/* 

MP2/6-311G(d,p) 

CH2Br2 

298.3 

361.1 

-123.0 

1.4 

MP2/6-3 1 1  G(2d,2p) 

298.7 

(293.3)b 

360.2 

-124.3 

1.2 

a  From  ref.  34 


b  From  ref.  40 


Table  IX  :  Summary  of  Arrhenius  parameters  calculated  at  different  levels  of  theory  including  the  tunneling  correction 
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Figure  1  :  MP2/6-311G(2d,2p)  optimized  geometry  for  the  transition  state  involved  in  the 
OH  +  CH2Br2  reaction 


Figure  2  :  Arrhenius  plots  for  OH  +  CH2Br2  ->•  H20  +  CHBr2  for  the  different  levels  of 
theory  used  in  this  study.  The  notations  mp2/2d2p  to  mp4/++3df3pd  in  the  legend  mean 
PMP2/6-311G(2d,2p)  to  PMP4SDTQ/6-311++G(3df,3pd). 


